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Weak lensing predictions for coupled dark energy 
cosmologies at non-linear scales 
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ABSTRACT 

We present non-linear weak lensing predictions for coupled dark energy models using 
the CoDECS simulations. We calculate the shear correlation function and error co- 
variance expected for these models, for forthcoming ground-based (such as DES) and 
space-based (Euclid) weak lensing surveys. We obtain predictions for the discrimina- 
tory power of a ground-based survey similar to DES and a space-based survey such 
as Euclid in distinguishing between ACDM and coupled dark energy models; we show 
that using the non-linear lensing signal we could discriminate between ACDM and 
exponential constant coupling models with /3q > 0.1 at 4a confidence level with a 
DES-like survey, and /3q ^ 0.05 at 5a confidence level with Euclid. We also demon- 
strate that estimating the coupled dark energy models' non-linear power spectrum, 
using the ACDM Halofit fitting formula, results in biases in the shear correlation 
function that exceed the survey errors. 
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1 INTRODUCTION 

After more than a decade of continuous improvements in 
the accuracy of cosmological observations - which has led 
to the establishment of a broadly accepted representation of 
our Universe known as the Concordance Cosmological Model 
(CCM) - we are now entering the epoch of precision cosmol- 
ogy. The great wealth of high-precision cosmological data 
expected throughout the next few years offers the exciting 
prospect of tightly constraining the parameters of the CCM 
or possibly detecting deviations from its predictions. 

The unprece dented resolution of the Planck satellite 
|Ade et al.l l201ll ) in measuring the angular temperature 
fluctuations and the polarization of the Cosmic Microwave 
Background (CMB) radiation will extract a significant 
amount of new information beyond that yielded by previ- 
ous CMB experiments, affording very tight constraints on 
the initial conditions of the Universe. 

At the other end of the cosmic expansion history, i, e. at 
low redshifts, present and future surveys measuring the clus- 
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tering of s ources dTeam et al.1 201 ll) and gravit ational lens- 
ing effects |Massev et al.ll2007l ; lTeam et al]201ll ) will greatly 



improve our knowledge of the spatial distribution of bary- 
onic and cold dark matter (CDM) in the local Universe, and 
of its time evolution, providing new tests of the gravitational 
instability processes driving the growth of cosmic structures. 

In particular, one of the most mysterious phenom- 
ena characterizing the low-redshift Universe is the appear- 
ance of a Dark Energy (DE) component capable of driving 
the observed present acceleration of the co s mic expansion 
i Riess et ail Il998l; iPerlmutter et all Il999l; ISchmidt et all 



119981 ; I Astier et alj|2006l ; iKowalski et alj|2008h . and required 
to also explain a number of other observations, e.g. the 
angular power spectrum of CMB temperature anisotropics 
l|Komatsu et alj 120091 , 1201 ll ), the evolution of the number 
counts of mass i ve galaxy cluste r s as a function of redshif t 
l|Borganil |2006| ; IVikhlinin et~"aH 120091 ; iMantz et ail l201Ch , 
the angular correlation of galaxies in large galaxy sur - 
veys (|Percival et al.l200ll ; [c"ole et alj20o5lReid et alj201Ch . 
or the observed scale of th e Baryon Acoustic Oscillations 
(BAO) |Percival et al.ll2010t ). In the standard ACDM model, 
such a DE component is identified with a cosmological con- 
stant A, a quantity with negative pressure and constant en- 
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ergy density throughout the whole expansion history of the 
Universe, and with no spatial fluctuations. 

This simple picture is very successful in reproducing a 
wide range of observational data. However, the cosmological 
constant scenario suffers from serious conceptual problems 
concerning the extremely small value of the constant DE 
density as compared to the typical densities of the early Uni- 
verse , known as the "fine tuning problem" (see e.g. lWeinbergl 
1989), and the apparent coincidence that it dominates 
over CDM only at relatively rece nt cosmological epochs , 
the "coincidence problem" (see e.g. iHuev fc Wandeltl l2006). 
In order to overcome these problems, alternative models 
based on the dyna mic evolution of a classical scalar field 
have been proposed ilWetterichlll988l ; iRatra fc Peeblesll 19881 ; 
lArmendariz-Picon et al.ll2000l ). Abandoning the simple pic- 
ture of a cosmological constant, however, necessarily requires 
us to consider and to include in our models of the Universe 
the presence of spatial fluctuations and of possible interac- 
tions of the new physical degree of freedom represented by 
the DE scalar field. 

It is in this context that models of interacting DE 
have been proposed as a natural extensi on of the mini- 
mally coupled dynamic scalar field scena rio l|Wetterichll 19951 ; 
lAmendolal l200ol ; H arrar fc RosenI 120071 ). Although an in- 
teraction of the DE scalar field wit h baryonic particles i s 
tightly constrained by observations (|Hagiwara et alj 12002). 
the same bounds do not apply to the case of a selective 
interaction be t ween DE and CDM, as first speculated by 
IDamour et al] l|l990l ). which has therefore received sub- 
stantial attention as a realistic competitor to the standard 
ACDM model. 

Various different forms of interactions between 
DE and CDM particles (including massive neutrinos) 
have been proposed and investigated in the litera- 



ture (as e.q. by Amendola 12004; 


Caldera-Cabral et alj 


20091; Pettorino & Baccigalurjil \20o4. 


Amendola et al.ll2008l; 


Boehmer et al.l 20101; Kovama et al. 


20091; Honorez et all 


2010|). and their impact on the linear growth of 



densi t y perturbations (see e.g. |Pi Porto fc Amendola! 
20081; ICaldera-Cabral et al.1 120091; IValiviita et al.1 |2008|; 



Maierotto et al.l |2010| ; IValiviita et all |2010| ; IClemson et al.1 
20111 ) and on the nonlinear regime of structure forma- 



tion ilMaccio et all |2003; IBaldi et all I2OI0I; iBaldil 1 2011 



Li fc Barrow 
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.201 lal ; IBaldi fc Pettorinol l20ld ; 111 I2OIC ; 
l2011bh has been extensively studied in re- 
cent years. For many such models, robust and realistic 
observational constraints on the interaction st rength have 
been derived based on CMB and LSS data (|Bean et al.l 



l2008l;lLa Vacca et al.ll2009l; IXiall2009h. local dynamical tests 
I Kesden fc Kamionk owski 20 061; iKeselman et al.ll2009h . and 



Lyman- a observables i|Baldi fc Vielll2010l ). Although these 
observational bounds have strongly restricted the allowed 
parameter space for interacting DE cosmologies, none of 
them has yet been able to rule out the model, or to un- 
ambiguously detect the presence of a DE-CDM interaction 
with compelling statistical significance. 

In this respect, exciting times are ahead of us, with the 
realistic possibility of exploiting the joint power of forthcom- 
ing high-precision cosmological observations to break many 
of the existing degeneracies between competing cosmologi- 
cal models and finally disentangle the distinctive features of 
alternative scenarios. Dark energy interactions will be one 



of the issues that can be tested, and so the next generation 
of cosmological data will possibly provide a real indication 
of the nature of the DE phenomenon. 

In the present paper, we examine the usefulness of weak 
gravitational lensing for discriminating between interacting 
dark energy models. We wish to show how the lensing sig- 
nal depends on the dark energy interaction, and whether 
this dependence is sufficiently strong that it could be de- 
tected with forthcoming lensing surveys. In particular, we 
will provide forecasts for the capability of future large Weak 
Lensing (WL) surveys -both a ground-based survey similar 
to the Dark Energy Survey (DESfl an d a space-based sur- 
vey, i.e. EUCLIE0 - to detect a DE-CDM interaction. Our 
particular focus in this paper is the non-linear regime, as 
this regime provides much of the power for lensing. To this 
end, we exploit the full non-linear matte r power spec trum 
as predicted by the CoDECS simulations (Baldi 2011bl ). the 
largest suite of self-consistent and high-resolution N-body 
simulations for interacting DE cosmologies to date. 

The paper is organized as follows: in Sec. [5] we describe 
the main features of the interacting DE models under in- 
vestigation; in Sec. [3] we discuss gravitational lensing in the 
context of interacting DE models, and in Sec. U we describe 
the methods used to compute the necessary nonlinear power 
spectra. In Sec. [S] we present the results of our analysis, giv- 
ing forecasts for forthcoming lensing surveys; we draw our 
conclusions in Sec. [6] 



2 COUPLED DARK ENERGY MODELS 

Coupled DE (cDE) models have been widely investigated in 
the literature concerning their cosmological background evo- 
lution as well as the behaviour of linear an d nonlinear den- 
sity perturbations in these models (see e.g. Amendola 2000, 
2004 iPettorino fc Baccigalupill2008l : |Pi Porto fc Amendolal 



20081 ; IBaldi et alj|20ld : iLi fc BarrowH2011al ; lBaldill201ll . and 
references therein). Here we only briefly introduce the defini- 
tions and the notation that will be assumed throughout the 
paper for the different cDE models; we refer the interested 
reader to the literature above for a thorough discussion of 
cDE scenarios. 

In the present work, we will consider cDE models de- 
fined by the following set of background dynamic equations: 



p b + 3Hp b = , 

pr + 4:Hpr = , 

2 1 

SH = T7T (Pr +Pc + Pb + , 



(1) 

(2) 

(3) 
(4) 

(5) 



where an overdot represents a derivative with respect to the 
cosmic time t, H = a/a is the Hubble function, V((f>) is 
the scalar field self-interaction potential, Mpi = l/V87rG is 
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the reduced Planck Mass, and the subscripts b ,c,r, indicate 
baryons, CDM, and radiation, respectively. 

The function /3 C (4>) sets the direction and the strength 
of the energy-momentum flow between the DE scalar field 
cf> and the CDM fluid, while the function V((j>) determines 
the dynamical evolution of the DE density. In the present 
work we will consider two possible ch oices for each of these 
two functions, namely an exponential ijLucchin fc Matarresd 
1 19851 : IWetterichll 1988 1 and a SUGRA prax fc Martinlll999l ) 
potential, 



EXP 
SUGRA 



V(4>) = Ae~ 
V{<t>) = 



(6) 
(7) 



where a is a positive constant and where for simplicity the 
field (j> has been expressed in units of the reduced Planck 
mass Mp\, as well as both a constant and an exponentially 
growing coupling function /3 C (</>): 



(8) 



characterized by f3\ = and /3i > 0, respectively. The 
most relevant difference between the exponential potential 
and the SUGRA potential relies on the fact that the 
latter features a global minimum at finite scalar field 
values; this allows for a change of direction of the scalar 
field motion, which is the main feat ure of the r ecently 
proposed "Bouncing cDE" scenario paldil 1201181 ). One 
should also notice that the notation introduced in Eqs. |l} 
51) corresponds to the original convention proposed by 



Amendolal (|2000h and has been adopted by several other 



studies, including the CoDECS project consider ed in the 
present work, but it differs by a constant factor v/2/3 from 
what is used in another part of the related literatu re (as 
e.g. iPettorino fc Baccigalupl 120081 ; iBaldi et alj|2010f ). The 
specific models considere d in t h e pres ent work h a ve bee n 
described in full detail bv lBaldil i|2011al ) and lBaldil (|2011bl ): 
we summarize them in Table [T] where the features and the 
specific parameters of each model are outlined. 

The evolution equations for linear density perturbations 
in the context of a cDE cosmology have been derived in the 
litera ture (see e.g. lAmen dola 2004; Pct torino fc Baccigalupil 
2008), and in the Newtonian limit of sub-horizon scales can 
be expressed as follows: 



S c = 



-2H 



l-Pc 



8 b = -2HS b + 4ttG [p b S b + p c 8. 



S c + 4nG [p b S b + p c 5 c T c ] , (9) 
(10) 



where S c ,b = 5p c ,b./pc,b are the relative density pertur- 
bations of the coupled CDM and uncoupled baryonic flu- 
ids, respectively, and where the scalar field dependence of 
the coupling function /3 c (<f>) has been omitted for simplic- 
ity. In Eq. (J9j), the factor T c = 1 + 4^(cf>)/3 represents 
an additional fifth-force mediated by the DE scalar field <j> 
for CDM perturbations, while the second term in the first 
square bracket at the right-hand-side of Eq. (£9|) is an ex- 
tra friction term on CDM fluctuations ari sing as a conse- 
quence of momentum conservation (see e.g. lAmendola| | 200l 
IPettorino fc BaccigalupHl2008l ; IBaldi et al.il2010l : iBaldaboil , 
for a derivation of Eqs. (I1I5I9I10P and for a detailed discus- 
sion of the extra friction and fifth force corrections to the 
evolution of linear perturbations). As a consequence of these 
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Figure 1. Power spectrum for ACDM and cDE models with con- 
stant coupling at z=0. 



two additional terms in the perturbed dynamic equations, 
CDM fluctuations will grow faster in cDE models with re- 
spect to a standard ACDM cosmology, thereby reaching a 
higher as normalization at z — if starting from the same 
amplitude at the last scattering surface zcmb ~ 1100, as 
shown in the last column of Table [T] However, in the non- 
linear regime the interplay between the friction term and 
the fifth force is not so straightforward as for the case of lin- 
ear perturbations, due to the fact that as a consequence of 
virialization processes, the local velocity field will not nec- 
essarily be aligned to the local gradient of the gravitational 
potential, as one can see from the three-dimentional general- 
ization of Eq. (|9} to a system of point-like massive particles: 



V c = Pc{<t>)—ptic — V 

V6 



(11) 



where r Ctb is the physical distance of the target coupled 
particle from the other CDM and baryonic particles, re- 
spectively. The effect of the friction term in the nonlinear 
regime has been shown to induce a suppression of small- 
scale power in the cDE models with respect to the nonlinear 
power that would be inferred based on the large-sc ale as 
normalization in the context of a ACDM universe paldil 
l201ll: iBaldil l2011bh . Such suppression will have important 
consequences on the weak lensing constraints on cDE mod- 
els that we want to address in the present work. Therefore, 
although it is possible to estimate the full matter power in 
cDE scenarios by applying nonlinear corrections (calibrated 
on ACDM simulations) to the re- normalized linear powe r 
spectrum (as recently done e.g. by lAmendola et al.ll201l1 ). 
in order to reach high accuracy at scales relevant for present 
and future large lensing surveys it is necessary to rely on 
a fully nonlinear treatment of cDE scenarios via specific N- 
body simulations. A discussion on the comparison between 
these two approaches is presented in Section [5.31 

Figure [1] shows the power spectra for each of the con- 
stant coupling (y8i = 0) models normalised by WMAP7. 
The values of these couplings were chosen since cDE mod- 
els with /3o ^ 0.15 can fit the angular diameter distance to 
decoupling measured by WMAP7, so these are of particu- 
lar interest as they are consistent with current observations 
of the background, but may on the other hand affect the 
growth of structures. It can be seen that there is a 2-7% 
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Model 


Potential 


a 


00 


/3i 


w^(z = 0) 


<Tg(z = 0) 


ACDM 


V{4>) = A 








-1.0 


0.809 


EXP001 


V(4>) = Ae- a ^ 


0.08 


0.05 





-0.997 


0.825 


EXP002 


V(4>) = Ae-^ 


0.08 


0.1 





-0.995 


0.875 


EXP003 


V(4>) = Ae- a <l> 


0.08 


0.15 





-0.992 


0.967 


EXP008e3 


V(<t>) = Ae" 01 ^ 


0.08 


0.4 


3 


-0.982 


0.895 


SUGRA003 


V(<j>) = A^^e^l 2 


2.15 


-0.15 





-0.901 


0.806 



Table 1. Interacting dark energy models considered in this work. In addition to the concordance ACDM model, we consider the 
exponential potential with three interaction strengths; the exponential potential with a time- varying strength; and the SUGR A potential. 



difference in the 2 = power spectrum between ACDM and 
EXP001, the lowest of the couplings investigated here, and 
a 25-65% difference between ACDM and the highest of the 
couplings, EXP003. 

In this analysis we do not use the simulated power spec- 
trum directly but instead use the ratio between the ACDM 
and cDE power spectra to find the difference in the growth 
of modes for different couplings with the same initial condi- 
tions. Using this method reduces the error associated with 
the limited number of independent fc-modes that enter the 
computation of the power in each k bin to only the error on 
the ACDM power spectrum. 



3 LENSING IN COUPLED DARK ENERGY 
COSMOLOGIES 

Now we present the framework for calculating the gravita- 
tional lensing signal in the cDE scenario. The way that light 
is deflected along the path from its source to an observer 
is determined by the mass distribution and the geometry of 
the Universe. The deflections of light lead to distortions of 
the observed image of the source. The mapping between the 
original source shape and the observed image is given by 



A = 



1 



k — 71 
—72 



-72 

K. + 71 



(12) 



(Bartclma nn fc Schn eider 20Ql|) where the convergence, ft, 
causes an isotropic dilation and the shear, 7 = 71 + 172, 
changes the ellipticity. k is challenging to measure, as the 
original size of the source is unknown; equally 7 cannot be 
measured for a single source as the intrinsic ellipticity of 
the source is unknown. However if the shear of a large num- 
ber of sources is correlated, then the lensing signal can be 
measured as a correlation function (insofar as the intrinsic 
ellipticities are not themselves correlated; see discussion in 
section [S] below). Therefore we will be interested in the shear 
correlatio n function in order to quant ify our predictions, 
given by (jBartelmann fc Sch neider 20q3) 



dl^-P K (l)J (W), 



(13) 



where 6 is the angular distance between the correlated 
sources, I is the angular w avenumber and the lensing power 
spect rum P K is given by |Bacon et al.l 120051 ; iMassev et al.l 
12007ft 



with the weight functions 
W( X ) = 



dx'G(x') 



(15) 



d X W 1 (x)W2{ X )a^r n {a) z P s ( l -, X ) ,(14) 



where x is comoving distance, xh is the comoving distance to 
the horizon and G(x) is the normalised distribution of the 
sources in comoving distance, corresponding to a redshift 
distribution for the sources. We use two weight functions in 
Equation [14] since we are using tomographic weak lensing. 
Equation (|15[) is valid for flat cosmologies, which are all that 
are considered in this paper. Usually Eq. (j 141) is written with 
the assumption f2 m (a) = f2 m /a' j ; however the form above 
does not include such an assumption, as coupling CDM and 
DE means that f2 m has a different dependence on time, as 
shown in Eqs. |T]|4]). 

We have modified the COSMOS Co s moMC 
code llLesgourgues et"al ] 120071 : iLewis fc Bridlel [2002; 
IMassev et al.ll2007l ). which calculates the combined shear 
correlation function from the theoretical power spectrum 
prediction given by CosmoMC, to include cross-correlation 
of redshift bins and to calculate the predicted weak lensing 
signal directly from the cDE model power spectra, according 
to Eqs. (|13H15p . We will now use these results to estimate 
the discriminatory power from lensing between different 
coupled DE models. 



4 SIMULATIONS 

For our analysis we will rely on the public nonlinear power 
spectrum da ta computed from the CoDECS simulations 
(Baldi 2011b), the largest suite of cosmological N-body sim- 
ulations for cD E models to date, c arried out with the mod- 
ified version bv lBaldi etakl l|2010t ) o f the widely us ed Tree- 
PM parallel N-body code GADGET (|SpringeW2005T ). In par- 
ticular we will consider the H-CoDECS suite that includes hy- 
drodynamical simulations of all the cDE models summarized 
in Table [T] on relatively small scales. More specifically, the 
H-CoDECS runs follow the evolution of 512 3 CDM and 512 3 
gas particles in a cosmological comoving box of 80 Mpc//i 
a side, with a mass resolution at z — of m c — 2.39 x 10 8 
M Q //i for CDM and m b = 4.78 x 10 7 M Q /h for baryons, 
and a force resolution set by the gravitational softening 
e 9 = 3.5 kpc/fe. Adiabatic hydrodynamical forces on the gas 
particles are computed by means of the entropy conserv- 
ing formulation of Smoothed Particle Hydrodynamics (SPH, 
ISpringel fc HernquistJ I2OO2T ) and other radiative processes 
such as gas cooling, star formation, or feedback mechanisms 
are not included in the simulations. 



Weak lensing in coupled dark energy cosmologies 5 



Parameter 



Value 



H 70.3 km s" 1 Mpc" 1 

J1 C DM 0.226 
U UE 0.729 

AOs) 2.42 x 10~ 9 (0.801 for ACDM) 

fl b 0.0451 

n s 0.966 



Table 2. The set of cosmological parameters at z = assumed 
for all the models included in the CoDECS project, consistent 
with t he latest results of the WMAP collaboration for CMB data 
alone llKomatsu et al.ll201lf) . 





n/ 


Area/ 


Survey 


galaxy arcmin 2 


degree^ 


DES 


13 


5000 


Euclid 


30 


15000 



Table 3. Galaxy density, n, and area assumed for our fiducial 
DES and Euclid surveys. 



Initial conditions are generated at z% = 99 by reseat- 
ing, with the appropriate growth factor for each specific 
model, the displacements obtained for a particular ran- 
dom field realization of the linear power spectrum Pn n (k) 
at zcmb- This power spectrum is com puted by the pub - 
licly available Boltzmann code CAMB (jLewis et alJ l2000l ) 
for a ACDM cosmology with parameters consistent with the 
latest "C MB only Maximum L ikelihood" constraints from 
WMAP7 (jKomatsu et al.1 120111 ). which are summarized in 
Table [2] This means that all the different simulations have 
exactly the same initial conditions at zqmb, and their differ- 
ent features at low redshifts depend uniquely on the different 
cosmology in place between last scattering and the present 
time. 

The H-CoDECS matter power spectra have been com- 
puted by evaluating the density of the different matter com- 
ponents on a grid with the same size of the PM grid (i. e. 512 3 
grid nodes) through a Cloud-in-Cell mass assignment of the 
different matter species and of the total matter distribution. 
This procedure allows us to compute the power spectrum 
up to scales corresponding to the Nyquist frequency of the 
grid, i.e. &N y = nN/L « 20.0/i/Mpc. Beyond this limit- 
ing frequency, the pow er spectrum has been computed with 
the fo lding method of I Jenkins et alJ (| 19981 ): IColombi et al.l 
(2008), and the two estimations have been smoothly inter- 
polated around fcNy Finally, the combined power spectrum 
has been truncated at scales where the shot-noise reaches 
10% of the measured power. 

With the power spectra computed with the procedure 
just described, we have investigated how future weak lensing 
probes could perform in constraining cDE cosmologies, as 
discussed in the next Section. 



5 RESULTS 

We calculated the combined shear correlation function for 
each of our models using equations 11341 151 We consider two 







DES 


Euclid 


Model 




Ax 2 


Ax 2 


EXP001 


0.05 


3 


30 


EXP002 


0.1 


18 


480 


EXP003 


0.15 


340 


3300 



Table 4. Best fit Ax 2 for different couplings, using errors calcu- 
lated for DES and Euclid surveys. 



types of survey: a ground-based survey modelled on DES, 
and a space-based survey, Euclid; the adopted galaxy density 
and survey area are shown in Table [3] In calculating the 
shear correlation function for these surveys we therefore use 
a DES-like redshift distribution given by 



n(z) 



(z a + z ab )/(z b + c) 



(16) 



where a — 0.612, b = 8.125, c = 0.62, and a space survey 
redshift distribution for Euclid given by 



n(z) 



cvEo— exp(- 
2, 



-(z/z ) 



(17) 



27 as used in 



where a = 2, = 3 /2, z = 0.63 and E 
iBevnon et ail (|2010T ). We also calculated simulated covari- 
ance matrices including sample va riance and shape nois e in 
a similar way to that calcula ted in IBevnon et all (|20ld ) us- 
ing the Horizon simulation (|Tevssier et al.l 20091 ): here we 
used 81 patches of 3.4 square degrees to estimate cosmolog- 
ical sample variance, and assumed an intrinsic shape noise 
of <t 7 = 0.2 on each component of the shear. 

In order to examine whether interacting dark energy 
models can be detected by forthcoming space and ground- 
based missions, we can assess the difference in x 2 between 
the best-fit ACDM and best-fit interacting DE model for a 
given dataset. One could choose a fiducial ACDM shear cor- 
relation function with realistic error-bars, and find the best- 
fit interacting DE model for this; but it is more convenient 
computationally to choose a fiducial interacting DE model 
and vary parameters of the easily obtained ACDM models 
to find the best standard cosmology fit. The difference in 
X 2 between the two best-fit models is the same whichever 
way round we choose, and is a measure of our ability to 
distinguish between the two types of model. 

We ran CosmoMC to find the best fit ACDM models for 
each of the cDE models with different CDM couplings. We 
used the following parameter space: fi m 0.5, 0.5 ^ 
a 8 «S 1, 0.4 «J h < 1, -2 «; w «; and 0.01 < Sl b < 0.15. 
The tomographic lensing results were studied for 3 cross- 
correlated redshift bins of equal size between z — 0.3 and 
z = 1.5 and 1' < 9 < 90'. 



5.1 Constant coupling models with an 
exponential potential 

In this section we look at how introducing a constant cou- 
pling between DM and DE (models EXP001-3 in Table [TJ 
affects the weak lensing signal. The shear correlation func- 
tions, with WMAP7 initial conditions, are shown in Figure 
[2] Note that /?o primarily changes the amplitude of the corre- 
lation function, with an additional slight alteration in slope. 
The difference in x 2 f° r each of the different constant cou- 
plings is shown in Table [4] and we see that lensing with 
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Figure 2. Correlation function predicted for cDE models with error estimates for DES (left) and Euclid (right) surveys using WMAP7 
best fit parameters. 



Euclid should be able to discriminate between /Jo ^ 0.05 
and ACDM at a confidence level of 5a, while DES should 
be able to discriminate between /Jo 0.1 and ACDM at a 
confidence level of 4rj. 

Figure [3] shows that the best fit ACDM models for each 
of the couplings occupy quite different parameter regions, 
especially for Euclid. The discrepancies between DES and 
Euclid predictions in these plots are found to be due to the 
off-diagonal covariance matrix terms; this can be seen by 
examining the best fit models for DES and Euclid along with 
the cDE model we are trying to fit. The best fit for our DES 
survey appears to be a worse fit at small 6 and a better fit at 
large 6 than the Euclid best fit. This is due to the covariance 
being largest for large angles and high redshifts. So while 
DES has a larger contribution from shape noise at small 8 
allowing a worse fit on small scales, conversely Euclid is more 
sensitive to the covariance on large scales. This descrepancy 
between the DES and Euclid best fit ACDM increases as /Jo 
increases. 

These results show that if dark energy and dark mat- 
ter truly do interact in the way described by our class of 
models, and we attempt to fit a ACDM cosmology to the 
observations, then we will infer increased values of Ho and 
as, and a decrease in w and fl m as /3p increas es. 

It shou l d be noted that iKirk et all (|201l] ) and 
iLaszlo et alj l|201lf) have recently shown that the effects of 
modified gravity and alternative dark energy models can 
be degenerate with systematics due to intrinsic alignments. 
Baryonic physics has also been shown to have possibly large 
effects on the mat ter power spectrum from scales as small 
as fc= 0.3 h/Mpc (|van Daalen et all l201ll ; ISemboloni et al.l 
l201ll ). In this paper we do not include these effects, as we 
are seeking to present the pure shear signal predictions. Our 
results should therefore be considered best-case predictions 
which will be diluted by the impact of systematic and bary- 
onic effects. 



5.2 Other potentials and coupling 

Although for the previous section we restricted ourselves to 
looking at constant coupling models with an exponential po- 
tential, the cDE model has the freedom to examine different 
potentials and an evolving coupling. Two of the CoDECS 





DES 


Euclid 


Model 


Ax 2 


Ax 2 


EXP008e3 


64 


570 


SUGRA003 


16 


100 



Table 6. Best fit Ax 2 for EXP008e3 and SUGRA003 using errors 
calculated for DES and Euclid. 

simulations explore this freedom: EXP008e3, which has the 
same potential as the models in the previous section but with 
an evolving coupling, and SUGRA003, which has a SUGRA 
potential with a constant coupling. Since there is not yet a 
suite of these types of simulations exploring the full range 
of parameter space, we have included them as lone exam- 
ples simply to demonstrate the range of the cDE model. 
The power spectrum for these models is shown in Figure 3] 
where we can see that for the EXP008e3 model we get simi- 
lar differences between the cDE model and ACDM to those 
shown in the larger constant coupling models (EXP002/3). 
On the other hand, the SUGRA003 model has smaller dif- 
ferences to this at large scales and much larger differences at 
small scales (almost 100% at k = Wh/Mpc) demonstrating 
how important it is to carry out full simulations of these 
models in order to obtain small scale predictions. 

We again attempted to find a best fit ACDM model us- 
ing CosmoMC and the \ f° r the best fit result, shown in 
Table [6j demonstrates that for these particular models we 
would be able to exclude both models at > la for both DES 
and Euclid. Further investigation of these types of model 
would allow constraints to be made on the parameters char- 
acterizing the coupling and the potential. 

5.3 Comparison of simulations and Halofit 

In section [2] we discussed the importance of using N-body 
simulations ov er using ACDM n on-linear fitting formulas 
such as Halofit (|Smith et al. 2003) to estimate the non-linear 
power spectrum for cDE models. In Figure [5] we show that 
the use of Halofit to estimate the non-linear power spectrum 
results in errors in the shear correlation function that exceed 
the statistical errors, for each of the surveys and for all of 
the models considered. This demonstrates the importance 
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Figure 3. Constraints on Q m , cs, n s , w and Ho- The light grey contours show the 68% and 95% confidence limits for DES, while the 
dark grey contours show the 68% and 95% confidence limits for Euclid. 



Survey 


Model 


w 


H 


CT8 




n 3 


DES 


EXPOOl 

EXP002 

EXP003 

EXP008e3 

SUGRA003 


-0.974 ±0.020 
-1.012 ±0.047 
-1.110 ±0.045 
-0.981 ±0.048 
-0.755 ±0.044 


69.2 ±3.5 
82.7 ±9.9 
95.1 ±2.8 

77.3 ± 10.0 
81.1±6.1 


0.834 ±0.005 
0.881 ±0.010 
0.946 ± 0.008 
0.889 ±0.010 
0.760 ±0.013 


0.264 ± 0.003 
0.259 ±0.005 
0.258 ± 0.004 
0.262 ±0.005 
0.305 ±0.008 


0.952 ±0.013 
0.973 ± 0.012 
0.947 ± 0.009 
0.954 ±0.014 
0.760 ± 0.013 


Euclid 


EXPOOl 

EXP002 

EXP003 

EXP008e3 

SUGRA003 


-0.974 ±0.020 
-0.888 ±0.020 
-1.004 ±0.020 
-0.881 ±0.020 
-0.804 ±0.020 


69.2 ±3.5 
66.1 ± 1.5 

73.3 ±1.3 
65.6 ±0.5 

85.4 ±2.2 


0.834 ±0.005 
0.918 ±0.004 
1.060 ±0.002 
0.935 ±0.004 
0.745 ± 0.004 


0.264 ± 0.003 
0.251 ±0.002 
0.218 ±0.001 
0.247 ±0.002 
0.314 ± 0.004 


0.952 ±0.013 
0.956 ± 0.018 
1.009 ± 0.007 
0.922 ± 0.016 
1.092 ± 0.007 



Table 5. Marginalised parameters for ACDM fit to models for DES and Euclid surveys with la errors. 
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Figure 5. Difference between shear correlation function calculated using simulations and shear correlation function calculated using 
Halofit. Also shown is the measurement error (from sample variance and shape noise) for DES (left) and Euclid (right) using WMAP7 
best fit parameters. 
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Figure 4. Power spectrum for an evolving coupling model with 
an exponential potential (EXP008c3), and a constant coupling 
model with a SUGRA potential. 

of using N-body simulations to predict the non-linear mat- 
ter power spectrum for cDE models, and that further sim- 
ulations for a variety of cDE models are required to make 
accurate weak lensing forecasts using non-linear scales. 



6 CONCLUSIONS 

In this paper we have presented weak lensing predictions for 
cDE models using the non-linear power spectrum calculated 
by the CoDECS simulations. 

We have calculated the total shear power spectrum for 
each of the models, and used CosmoMC to find the best 
fit ACDM model; we have demonstrated the discriminatory 
power of future lensing surveys such as DES and Euclid, 
where it should be possible to tightly constrain constant 
coupling models with exponential potentials to j3o < 0.05 
with Euclid, or /3o < 0.1 with DES. However, this should 
be considered a best-case scenario, since the inclusion of 
intrinsic alignments and baryonic physics may impact the 
constraining power; this will be the subject of future work. 

We have shown that for cDE models with larger cou- 
pling there is a clear difference between the best fit ACDM 
for the same model but different surveys. This difference is 



due to the dominance of the off-diagonal covariance matrix 
terms over the diagonal for larger surveys, and shows the 
importance of including these off-diagonal terms in weak 
lensing predictions. 

We have also calculated the expected signal for a non- 
constant coupling model and a non-exponential potential 
model. These models could be excluded by 2a for a DES- 
like survey and > 7 a for Euclid. However we have not ob- 
tained constraints on the parameters of these types of model, 
since currently N-body simulations for these models have 
only been run with one parameter set. A substantial set of 
simulations would be required in order to properly sample 
the parameter space of these more complex scenarios. This 
will be a worthwhile task, as the effects of these cosmolo- 
gies appear to be more difficult to detect in the background 
and in the linear regime with respect to standard interacting 
dark energy models, making non-linear N-body simulations 
vital for realistic lensing predictions. 

We have also shown the size of the error on weak lens- 
ing predictions if a ACDM non-linear fitting formula, such 
as Halofit, is used to estimate the matter power spectrum, 
instead of using simulations. We find that this Halofit er- 
ror is larger than the statistical error for the DES and Eu- 
clid surveys, and for all the models considered here. This 
demonstrates the importance of using a full N-body code to 
estimate the non-linear power spectrum. 



ACKNOWLEDGEMENTS 

DB is supported by an RCUK Academic Fellowship. KK 
is supported by the STFC (grant no. ST/H002774/1), a 
European Research Council Starting Grant and the Lev- 
crhulme trust. EB is funded by an STFC PhD studentship. 
MB acknowledges support by the DFG Cluster of Excellence 
"Origin and Structure of the Universe" and by the TRR33 
Transregio Collaborative Research Network on the "Dark 
Universe" . 

REFERENCES 

Ade P. A. R., et al., 2011 



Weak lensing in coupled dark energy cosmologies 9 



Amendola L., 2000, Phys. Rev., D62, 043511 
Amendola L., 2004, Phys. Rev., D69, 103524 
Amendola L., Baldi M., Wetterich C, 2008, Phys. Rev., 
D78, 023015 

Amendola L., Pettorino V., Quercellini C, Vollmer A., 
2011 

Armendariz-Picon C, Mukhanov V. F., Steinhardt P. J., 

2000, Phys. Rev. Lett., 85, 4438 
Astier P., et al, 2006, Astron. Astrophys., 447, 31 
Bacon D. J., et al., 2005, Mon. Not. Roy. Astron. Soc., 363, 

723 

Baldi M., 2011, Mon. Not. Roy. Astron. Soc., 414, 116 
Baldi M., 2011a, Mon. Not. Roy. Astron. Soc. submitted 

[arXiv: 1107.5049] 
Baldi M., 2011b. larXi"v:1109.5695l 

Baldi M., 2011, Mon. Not. Roy. Astron. Soc, 411, 1077 

Baldi M., Pettorino V., 2010. larXiv:1006.3761l 

Baldi M., Pettorino V., Robbers G., Springel V., 2010, MN- 

RAS, 403, 1684 

Baldi M., Viel M., 2010. larXiv:1007.3736l 

Bartelmann M., Schneider P., 2001, Phys. Rept., 340, 291 

Bean R., Flanagan E. E., Laszlo I., Trodden M., 2008, Phys. 

Rev., D78, 123514 
Beynon E., Bacon D. J., Koyama K., 2010, Mon. Not. Roy. 

Astron. Soc, 403, 353 
Boehmer C. G., Caldera-Cabral G., Chan N., Lazkoz R., 

Maartens R., 2010, Phys. Rev., D81, 083003 
Borgani S., 2006 

Brax P., Martin J., 1999, Phys. Lett., B468, 40 
Caldera-Cabral G., Maartens R., Schaefer B. M., 2009, 

JCAP, 0907, 027 
Caldera-Cabral G., Maartens R., Urena-Lopez L. A., 2009, 

Phys. Rev., D79, 063518 
Clemson T., Koyama K., Zhao G.-B., Maartens R., Valivi- 

ita J., 2011 

Cole S., et al., 2005, Mon. Not. Roy. Astron. Soc, 362, 505 
Colombi S., Jaffe A. H., Novikov D., Pichon C, 2008, 
larXiv:0811.0313l 

Damour T., Gibbons G. W., Gundlach C, 1990, Phys. Rev. 
Lett., 64, 123 

Di Porto C, Amendola L., 2008, Phys. Rev., D77, 083508 
Farrar G. R., Rosen R. A., 2007, Physical Review Letters, 
98, 171302 

Hagiwara K., et al., 2002, Phys. Rev., D66, 010001 
Honorez L. L., Reid B. A., Mena O., Verde L., Jimenez R., 

2010, 1 arXiv: 1006.08771 
Huey G., Wandelt B. D., 2006, Phys. Rev., D74, 023519 
Jenkins A., et al., 1998, Astrophys. J., 499, 20 
Kesden M., Kamionkowski M., 2006, Phys. Rev., D74, 

083007 

Keselman J. A., Nusser A., Peebles P. J. E., 2009, Phys. 

Rev., D80, 063517 
Kirk D., Laszlo I., Bridle S., Bean R., 2011 
Komatsu E., et al., 2009, Astrophys. J. Suppl., 180, 330 
Komatsu E., et al., 2011, Astrophys. J. Suppl., 192, 18 
Kowalski M., et al., 2008, Astrophys. J., 686, 749 
Koyama K., Maartens R., Song Y.-S., 2009, JCAP, 0910, 

017 

La Vacca G., Kristiansen J. R., Colombo L. P. L., Mainini 

R., Bonometto S. A., 2009, JCAP, 0904, 007 
Laszlo I., Bean R., Kirk D., Bridle S., 2011 



Lesgourgues J., Viel M., Haehnelt M. G., Massey R., 2007, 

JCAP, 0711, 008 
Lewis A., Bridle S., 2002, Phys. Rev., D66, 103511 
Lewis A., Challinor A., Lasenby A., 2000, Astrophys. J., 

538, 473 
Li B., 2010 

Li B., Barrow J. D., 2011a, Phys. Rev., D83, 024007 
Li B., Barrow J. D., 2011b, Mon. Not. Roy. Astron. Soc, 
413, 262 

Lucchin F., Matarrese S., 1985, Phys. Rev., D32, 1316 
Maccio A. V., Quercellini C, Mainini R., Amendola L., 

Bonometto S. A., 2004, Phys. Rev., D69, 123516 
Majerotto E., Valiviita J., Maartens R., 2010, Mon. Not. 

Roy. Astron. Soc, 402, 2344 
Mantz A., Allen S. W., Rapetti D., Ebeling H., 2010, Mon. 

Not. Roy. Astron. Soc, 406, 1759 
Massey R., et al., 2007, Astrophys. J. Suppl., 172, 239 
Percival W. J., et al., 2001, Mon. Not. Roy. Astron. Soc, 

327, 1297 

Percival W. J., et al., 2010, Mon. Not. Roy. Astron. Soc, 
401, 2148 

Perlmutter S., et al., 1999, Astrophys. J., 517, 565 
Pettorino V., Baccigalupi C, 2008, Phys. Rev., D77, 
103003 

Ratra B., Peebles P. J. E., 1988, Phys. Rev., D37, 3406 
Reid B. A., et al., 2010, Mon. Not. Roy. Astron. Soc, 404, 
60 

Riess A. G., et al., 1998, Astron. J., 116, 1009 
Schmidt B. P., et al., 1998, Astrophys.J., 507, 46 
Semboloni E., Hoekstra H., Schaye J., van Daalen M. P., 

McCarthy I. J., 2011. larXiv:1105. 10751 
Smith R. E., et al., 2003, Mon. Not. Roy. Astron. Soc, 341, 

1311 

Springel V., 2005, Mon. Not. Roy. Astron. Soc, 364, 1105 
Springel V., Hernquist L., 2002, Mon. Not. Roy. Astron. 

Soc, 333, 649 
Team E., et al., 2011 

Teyssier R., et al., 2009, Astron. Astrophys., 497, 335 
Valiviita J., Maartens R., Majerotto E., 2010, Mon. Not. 

Roy. Astron. Soc, 402, 2355 
Valiviita J., Majerotto E., Maartens R., 2008, JCAP, 0807, 

020 

van Daalen M. P., Schaye J., Booth C, Vecchia C. D., 2011 
Vikhlinin A., et al., 2009, Astrophys. J., 692, 1060 
Weinberg S., 1989, Rev. Mod. Phys., 61, 1 
Wetterich C, 1988, Nucl. Phys., B302, 668 
Wetterich C, 1995, Astron. Astrophys., 301, 321 
Xia J.-Q., 2009, Phys. Rev., D80, 103514 



